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Abstract 

The behavior of orbits for iterated logistic maps has been widely studied since the dawn of 
discrete dynamics as a research field, in particular in the context of the complex family /: C —)■ C, 
parametrized as fdz) = + c, with c G C. However, little is is known about orbit behavior if 

the map changes along with the iterations. We investigate in which ways the traditional theory 
of Fatou-Julia may still apply in this case, illustrating how the iteration pattern (symbolic 
template) can affect the topology of the Julia and Mandelbrot sets. 

We briefly discuss the potential of this extension towards a variety of applications in genetic 
and neural coding, since it investigates how an occasional or a reoccurring error in a replication 
or learning algorithm may affect the dynamic outcome. 


1 Introduction 

1.1 Discrete dynamics of the logistic family 

The family of logistic maps has been over the years one of the most studied examples in the theory 
of discrete dynamical systems. In the context of real interval maps, typically parametrized as 
/^: [0,1] —>■ [0,1], f^{x) = fix{l — x), for fi G [0,4], results from kneading theory classify the 
possible orbits, and provide a relationship between the combinatorics of critical orbits and the 
complexity of the corresponding map (measured, for example, via its topological entropy). In the 
context of complex functions /c: C —)■ C, fc{z) = + c, for c G C, results going all the way back to 

the original theory of Fatou and Julia describe thoroughly the behavior of the orbits in the dynamic 
complex plane, as well as phenomena in the parameter plane. 

On the practical side, discrete iterations in general, and the dynamics of quadratic functions in 
particular, have been used to model natural processes. For example, the quadratic family has been 
used for more than a decade to model integrate and hre neurons [8, 3], and iterations of simple 
discrete maps are the ideal candidate for modeling replication in genetic algorithms. 

However, in nature, it is unlikely that systems evolve according to the same identical dynamics 
along time. Rather, one expects that occasional, or even periodic errors may be made in the iteration 
process, or even that the iteration scheme may change in time, according to the system’s new needs 
as it is adapting. Therefore, a more realistic mathematical context to model such phenomena 
is to consider time-dependant iterations, in which the iterated map may change between steps, 
evolve in time, or appear (with variable frequency) in conjunction with other maps in the iteration. 
Based on this premise, some recent studies have looked at how the traditional theory extends 
when alternating two different real [1, 14] or complex maps [5, 6] in the iteration. In the latter 
references, the authors introduce the extended Julia sets corresponding to the alternated iteration 
of two complex quadratic maps. They show, analytically and numerically, that the dichotomy in 
the classical single map case (where the Julia sets are either totally connected, for some values of 
c, or totally disconnected, for the remaining values of c) does not hold in this case, and that Julia 
sets for alternated maps can also be disconnected without being totally disconnected. They furhter 
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relate this extension to the Fatou-Julia theorem in the case of complex polynomials of degree greater 
than two, and show that alternated Julia sets exhibit graphical alternation. 

We further extend this idea by looking at iteration of two different functions, fcg and f^, 
according to a more general binary symbolic sequence, in which the zero positions correspond 
to iterating the function /eg and the one positions correspond to iterating the function fd- This 
generalization is a more appropriate representation of replication or learning algorithms that appear 
in nature, with patterns that evolve in time, and which may involve occasional, random or periodic 
“errors”. In this paper we investigate, primarily from a visual and numerical perspective, the 
questions that arise for such complex iterations, and to what extent existing results are expected 
to hold in this more general setting. We are currently working on a theoretical frame to support 
analytical proofs for the type of results conjectured here. 

1.2 The traditional Julia and Mandelbrot sets 

The prisoner set of a complex map / is defined as the set of all points in the dynamic plane, whose 
orbits are bounded. The escape set of a complex map is the set of all points whose orbits are 
unbounded. The Julia set of / is definded as their common boundary J(/). The filled Julia set is 
the union of prisoner points with their boundary J(/). 

For polynomial maps, it has been shown that the connectivity of a map’s Julia set is tightly 
related to the structure of its critical orbits (i.e., the orbits of the map’s critical points). Due 
to extensive work spanning almost one centry, from Julia [11] and Fatou [9] until recent develop¬ 
ments [2, 13], we now have the following result, known as the Fatou-Julia Theorem: 

Theorem 1.1. For a polynomial with at least one critical orbit unbounded, the Julia set is totally 
disconnected if and only if all the bounded critical orbits are aperiodic. 

It was later shown [4, 7] that, as a consequence for the degree two case, Julia sets are either 
connected (if the orbit of the critical point 0 is bounded) or totally disconnected (if the orbit of 0 
is unbounded). 

If one had hoped similar results to hold for periodic iterations of more than one map, Danca 
et al. [6] showed that this dichotomy is broken even at the level of a simple alternation of two 
quadratic maps: aside from the two classical situations, the alternation also produced disconnected 
Julia sets which were not totally disconnected. While for periodic iterations one may still be able to 
build upon existing work and relate the structure of Julia sets with the critical alternating orbits, 
for nonperiodic iterations, however, the foundation is unlikely to remain applicable, and one may 
have to construct entirely new methods. 

1.3 Our extension: definitions and notations 

As with the traditional Julia set, we will be working with the complex quadratic family 

{fc : C —>■ C / fc{z) = z'^ + c, with c G C} 

However, for each iteration process, we will be using a pair of maps in this family, /eg and fd, as 
follows: 

Definition 1.2. Fix co,ci G C, and a sequence s = (sn)n>o £ 1}) (which we will call the 

symbolic template of the iteration). For any £ C, the template orbit Os(^o) = (?n)n>o is the 
sequence constructed recursively, for every n > 0, as: 

?n+l = fcsn 
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In other words, for each n > 0, we iterate fcoiz) = z'^ + cq if the corresponding entry = 0 in the 
symbolic template, and we iterate fci{z) = + ci otherwise (i.e., if Sn = !)• The Julia set for a 

symbolic template system is naturally defined as an extension of the traditional one: 

Definition 1.3. Fix cq , ci G C. Then we define, for any s G >C({0,1}); 

The template prisoner set: Pco,ci(s) = {'^o G C / Os(.^o) 'Is bounded } 

The template escape set: £'co,ci(s) = {Co G C / Os(Co) 'not bounded } 

The template Julia set: Jcq,ci(s) = 9Pco,ci(s) = cl-Eco,ci(s) 

Since, in order to construct a template Julia set, we need both a complex parameter pair (co,ci) 
and a symbolic template, we can consider phenomena in two different parameter spaces: in (for 
a fixed template s) and in the template space of binary sequences /({0,1}) (for a fixed pair (cq, ci)). 
We define two types of “Mandelbrot” sets, as follows: 

Definition 1.4. Fix s G T({0,1}) symbolic sequence. The corresponding fixed-template Man¬ 
delbrot set is defined as: 


AJs = {(co,ci) G / Os(0) is bounded } 

Definition 1.5. Fix (co,ci) G C^. The corresponding fixed-maps Mandelbrot set is defined as: 

Alco.ci = (s G >C({0,1}) / Os(0) is bounded } 

Properties of the fixed template Mandelbrot set have been previously studied for period two tem¬ 
plates (i.e., alternating maps). In [5], the authors illustrate their corresponding extension of the 
Mandelbrot set as a subset of the quaternion filed, as well as some of its two-dimensional cross- 
sections. In the context of our extension, one may be interested to investigate a few directions, 
among which: 

- Assessing the topological (e.g., connectedness) and fractal properties (e.g., Hausdorff dimen¬ 
tion of the boundary) of hxed template Mandelbrot sets Alg) and how these vary with chang¬ 
ing the template. 

- Understanding topological and measure theoretical properties of fixed map Mandelbrot sets 
J^co.ci) and how there depend on changing the parameters cq and ci. 

~ Studying whether / which aspects of the Fatou-Julia theory still hold for this extension. For 
example, it is clearly no longer true that an unbounded template orbit of zero implies a totally 
disconnected template Julia set. However, is it possible that the Julia set remain connected 
if the template orbit of 0 is bounded? In which case , do we still have a result such as the 
following: 

Conjecture 1.6. For any template s, the fixed template Mandelbrot set Aig s equivalent with: 

A4s = {(co,ci) G / Jcg,ci(s) 'is connected } 

For any (co,ci) G C^, the fixed map Mandelbrot set AJco.ci is equivalent with: 

-^co,ci = (s G T({0,1}) / Jco,ci(s) is connected } 
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To the best of our knowledge, properties of the fixed-map Mandelbrot set have never been studied, 
and we believe this to be the first time the concept is being used. 

In the following section, we investigate iterations of complex qudratic maps according to a variety 
of templates, and we discuss in particular some questions that arise as natural extensions of the 
traditional problems. The paper is organized as follws: In Section 2 we discuss properties of our 
Julia and Mandelbrot sets for fixed periodic templates. We focus in particular on coupling an 
arbitrary quadratic map fc{z) = + c with the trivial map fo{z) = (whose Julia set is the unit 

circle). We observe how the connectivity patterns of the fixed template Julia and Mandelbrot sets 
change when the template itself is modified, to incorporate the two iterations in a different periodic 
mixture. In Section 3, we examine the same properties for nonperiodic templates. In Section 3.2 we 
view the insertion of the second map as an “error” in the iteration process of the first map, study 
the effects of such an error propagating along the template on the structure of the corresponding 
Julia set. In Section 4, we view the fixed map Mandelbrot set as a subset of [0,1]^, and we define 
hybrid Mandelbrot sets. Finally, in Section 5, we briefly discuss potential applications of this theory 
to the modeling of natural systems. 

2 Periodic template dynamics 

2.1 Periodic template Julia sets 

In this section, we will consider /c-periodic templates: s = [si,..., Sfc], with Sk+j = Sfc, for all j < 0. 
In previous work, Danca et al. have considered iterations of alternating quadratic maps [5, 6]. 
In our context, this corresponds to considering the two possible symbolic templates of period two 
s = [01] and s = [10]. 

In the references, the authors show that some basic properties and results are still inherited 
from the traditional single map case, but also that some new results emerge. For example, for any 
alternating orbit, the sequence of even and odd iterates are simultaneously bounded or unbounded, 
so that the Julia set of the alternated maps fcQ and fc^ and the Julia set of the quartic map o 
have the same connectivity type for any co,ci G C. The authors further investigate numerically, 
graphically and analytically the relationship between the connectivity type of Julia sets and the 
boundedness of the critical orbits of alternating maps, verifying the first part of the Fatou-Julia 
theorem in this case. On the other hand, they also show that the dichotomy (connected vs, totally 
disconnected) characteristic to classical Julia sets no longer applies under the extension. Other 
questions, such as conditions for local connectivity of the 4-dimensional Mandelbrot set or of 2- 
dimensional Mandelbrot slices, remain open to further investigation. 



Figure 1: Classical Julia sets for four different values of the parameter c, for which the sets are 
connected: A. c = —0.75; B. c = —0.117 — 0.76f; C. c = —0.62 — 0.432z; D. c = —0.117 — 0.856f. 
We will be using these values in our subsequent illustrations, for combinations of two such functions 
iterated along a symbolic template. 
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Similar properties are likely to extend to periodic template iterations for periods k > 2. Notice 
that a j-shift a on the periodic template cr(s) = cr[si,...Sfc] —)■ [si+j, ...Sfc+j] = translates as a 
polynomial transformation of degree j on the Julia set: 

One may look at classes defined by shifts on periodic templates, so that Julia sets within the same 
class are all related by such transformations. The connectedness properties of a Julia set may be 
preserved or not by the transformation, depending on the parameters cq and ci, and on the critical 
orbit properties produced by their combination. 

For templates of period 3 for example, there are 3 such classes (or types): 

(0.) All zero entries: s = [000] 

(1.) A single one entry: s = [001], s = [010], s = [100] 

(2.) Two one entries: s = [011], s = [101], s = [110] 

(3.) All one entries: s = [111] 

While higher period template Julia sets show similar symmetry properties as observed in the 
alternating case [5], new questions add to the discussion for higher periods. These refer to un¬ 
derstanding how the structure of the periodic template block affects the resulting Julia set, in 
combination with the parameter pair (co,ci). Connectivity of the Julia set is now influenced by 
three factors: (i) the maps cq and ci, (ii) the balance of how often one map is iterated versus the 
other (i.e., the number of Is vesus Os in the template block) and (hi) the location of these Is and 
Os along the block. 

Figures 2, 3 and 4 show, for example, template Julia sets for 3-periodic templates, for combinations 
of the map cq = 0 with the maps ci = —0.117 — 0.76i, ci = —0.62 — 0.432i ci = —0.5622 — 0.62i 
and respectively ci = —1 — 0.55i. For the first two of these (cq, ci) pairs, the classical Julia sets are 
shown in Figure 1; the other two Julia sets are totally disconnected. 

The template type seems to affect connectivity of the template Julia set as significantly as the 
parameters (co,ci). Clearly, it is not surprising that, for fixed cq and ci, using a different template 
type (e.g., 1 versus 2 ones in the 3-periodic case) will affect the connectivity of the Julia set. 
However, the results may be rather counterintuitive (see Figure 5). In addition to the template 
type, changing the position of the Is along the template may or may not affect connectivity, as 
shown in Figures 2 and 3 for the three different template blocks of type 2. In particular, the folding 
produced by applying a shift to the template may break or merge connectivity loci, depending on 
the template and on the complex parameters cq and ci used. 

For Cl = —0.62 — 0.432z (Figure 2), ci = —0.117 — 0.76z (Figure 3) and ci = —0.5622 — 0.62i 
(Figure 4) in combination with cq = 0, the Julia set becomes more disconnected as a result of 
increasing the contribution of ci versus cq. Indeed, in the first case, the template Julia set, connected 
for all cases of type one template (a single one and two zeros) remains connected (although a lot 
more complex) for all corresponding type 2 templates (two ones and one zero). In the second case, 
while all type 1 templates render connected Julia sets, these become disconnected for all cases of 
type 2 template. In the third case, the Julia set for type 1 templates, already disconnected, turns 
to dust for all type 2 templates. This is somewhat expectable, given the shape of the Julia set for 
each of these maps considered in isolation (the unit circle for c = 0, the Julia set shown in Figure Ic 
for c = —0.62 — 0.432i, in Figure lb, for c = —0.117 — 0.76z, and respectively a totally disconnected 
Julia set for c = —0.5622 — 0.62i). However, for ci = — 1 — 0.55i (Figure 5), the opposite happens: 
the template Julia sets are disconnected, but not totally disconnected for type 2 templates, and 
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Figure 2: Period 3 template Julia sets for cq = 0 and ci = —0.62 — 0.432z, for templates of 
type 1 (top row) and type 2 (bottom row). A. s = [Oil]; B. s = [101]; C. s = [110]; D. s = [001]; 
E. s = [010]; F. s = [100]. 



Figure 3: Period 3 template Julia sets for cq = 0 and ci = —0.117 — 0.76z, for templates of 
type 1 (top row) and type 2 (bottom row). A. s = [011]; B. s = [101]; C. s = [110]; D. s = [001]; 
E. s = [010]; F. s = [100]. 


become dust for type 1 templates, where the more substantial contribution of the map /co(^) = 
would suggest otherwise. 
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Figure 4: Period 3 template Julia sets for cq = 0 and ci = —0.5622 — 0.62i, for templates of 
type 1: A. s = [001]; B. s = [010]; C. s = [100]. The Julia sets for the same parameter values and 
templates of type 2, i.e. s = [Oil], s = [101], s = [110], are totally disconnected (dust). 



Figure 5: Period 3 template Julia sets for cq = 0 and ci = —1 — 0.55i, for templates of type 
2: A. s = [Oil]; B. s = [101]; C. s = [110]. The Julia sets for the same parameter values and 
templates of type 1, i.e. s = [001], s = [010], s = [100], are totally disconnected (dust). 

2.2 Periodic template Mandelbrot sets 

We next observe the structure of template Mandelbrot sets for periodic templates. We are inter¬ 
ested in particular to illustrate how the fractal structure of the set (represented by its Hausdorff 
dimension) and its connectedness properties change when changing the template density or regu¬ 
larity. 

Since, as in the original reference, we are using two iterated maps in combination, the template 
Mandelbort set can be seen as a 2-dimensional complex object (the locus of (cq, ci) G for which 
the template orbit of zero is bounded). One can visualize the connectedness of the MAndelbrot set 
for a fixed template by looking at a latice of 2-dimensional slices, each representing the behavior 
with respect to ci for a fixed cq (see in Figure 17 in Appendix 2). 

When we tracked the fractality of the boundary of template Mandelbrot sets, we found consis¬ 
tent trends, which we describe visually in Figure 6, for periodic templates, and in Figure 10, for 
nonperiodic templates, using consecutive zooms into various regions along the boundary, 
we show consecutive zooms into the boundary structure, at deeper and deeper levels). 
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Figure 6: Template Mandelbrot slice for the periodic template [Oil], for fixed cq = 
—0.2 + 0 .with two progressively zoomed in windows: one with low, the other with high Hausdorff 
dimension along the boundary. Three consecutive zoom-ins of the first instance are represented on 
top, and four zoom-ins of the second situation are represented on the bottom. The pattern persists 
at higher and higher levels, suggesting an alternation of high and low dimension intervals along the 
boundary of the set. 

3 Nonperiodic templates 

The case of periodic templates represents only a first extension, with basic properties expected 
to replicate quite naturally the theory for iterations or alternated maps. We next consider the 
more general case, of binary templates which are not necessarily periodic. In this case, the existing 
methods cannot be applied directly. To begin with, in this more general case, there is no extension 
for the classic concept of critical orbits, hence no corresponding possibility of encoding the behavior 
of the Julia set in the behavior of these orbits. As a consequence, results such as an extension of 
the Fatou-Julia theorem, or local connectivity of fixed-template Mandelbrot sets may be harder to 
determine, and require new methods. In this paper we only describe the setup, draw comparisons 
with the classical case and lay out a few of the interesting questions and conjectures that arise in 
this new framework. A more rigorous analysis, discussed briefly in the last section, is the subject 
of our current research. 

3.1 Julia sets for nonperiodic templates 

In Figure 7, one nonperiodic template was used to create the Julia set for three different combina¬ 
tions of parameters. Perhaps the first thing one notices is that the Julia sets for random templates 
still exhibit complex (fractal) structure, symmetry and alternations. However, we need to recall 
that, in all numerically generated figures, one can only use truncated representations of infinite 
templates, retaining only a specific number of iterations when representing the Julia set (for exam¬ 
ple, in all our figures we used 200 iterations). A natural question concerns the variability with the 
number of iterations (i.e., length of the truncated template) of the resulting approximation for the 
Julia set. 

In Figure 8, three different nonperiodic (randomly generated) templates were used to create the 
corresponding Julia sets, for the now familiar combination of parameter values ci = —0.62 — 0.432i 
and Co = 0. Notice that different templates introduce differences in the Julia set as significant as 
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Figure 7: Template Julia sets for one fixed nonperiodic template, and different pairs (cq, ci); A. 
Cl = —0.75 and cq = —1.2; B. ci = —1.2 and cq = —0.75; C. ci = —1.2 and cq = —0.8. For the 
simulation, the template ereating the Julia set was truncated to 200 iterations. 



Figure 8: Template Julia sets for three different randomly generated templates (from left to 
right), with parameters cq = 0 and, from top to bottom: ci = —0.75 (top), ci = —0.8 (middle) and 
Cl = 0.375 + 0.333i (bottom). For the simulation, the template creating the Julia set was truncated 
to 200 iterations. 

changes in the parameters (cq, ci). We then fixed Template 1 from Figure 8 and cq = 0, ci = — 0.8 
(that is, we started with the Julia set in the middle left panel). Using the same (co,ci) pair, we 
computed the Julia sets for a collection of templates of length 200 with identical entries up to the 
10th position. The variations among the Julia sets we obtained this way are clearly much smaller, 
as illustrated in Figure 9. To better understand this dependence, we consider the following: 
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Definition 3.1. For any n G N, we call the k-root of a template s = (sk) the finite sequence 
# = Sl, . . . ,Sfc. 

Definition 3.2. We say that two templates have the same k-root if they agree up to their k-th 
position. 

Based on our numerical simulations (see Figure 9), we propose the following conjectures: 

Conjecture 3.3. Fix a parameter pair (cq, ci), and a template s. For any <5 > 0, there exists n G N 
such that, if Sl and S 2 are two templates with a common n-root, then 


d{Jco,ci (si) 5 Jco,ci (S2)) < 5 

where d represents the Flausdorff distance between two sets. 

Conjecture 3.4. Fix a parameter pair (co,ci). For any 5 > 0, there exists N gN such that: 

d{Jco,ci{s^) - dco,ci(s)) <5 , for n>N 

In other words, the truncated Julia sets converge to the true Julia set in Hausdorff distance, as the 
root of the template is getting larger. 



Figure 9: Template Julia sets for ci = —0.8 and cq = 0, for three random truncated templates of 
length 200 with the same root of length / = 10. The first template (left) is the same as Template 1 
in Figure 8, hence the corresponding Julia set is identical with that in the middle left panel in the 
referenced figure. 


3.2 Propagating a perturbation 

One way to interpret the combination of two maps in the iteration scheme is to think of as the 
desired map, and of /cq as an “error,” or “perturbation” in the desired iteration. In this context, 
an all 1 template corresponds to a “perfect” replication process, in which a specific map is 
iterated identically for a (large or infinite) number of times, and any number of zero entries in the 
template correspond to as many intrusions of the erroneous map /cq in this iteration process. This 
intrusion can be periodic (as discussed in Section 2), or may occur at a random sequence of steps 
(as discussed in Section 3.1). In this section, we illustrate the changes in the Julia set produced by 
propagating a single error along a perfect template. 

In Appendix 1, we illustrate the effects on the Julia set of a single propagting error, for a few 
combinations of maps. In Figures 13 and 14, the desired map fc^ is combined with the trivial 
map fo{z) = , while in Figure 15, the error is a different (and unrelated) function, with complex 
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parameter cq 7 ^ 0. In Figure 16, the error is a small perturbation of the original function, i.e. 
Co = Cl + e, with a small e G C. 

In general, we notice that the impact of the same error on the Julia set is more substantial when 
the error occurs earlier in the iteration process. In fact, the Julia sets start to be reminiscent of the 
correct Julia set less than 50 iterations through the template length. Along the way, however, the 
connectivity and symmetry of the erroneous Julia changes at each step, sometimes in unexpected 
and counter-intuitive ways. Notice, for example, that the small perturbation to ci (in Figure 16) 
has at each step an effect on the Julia set which is comparable with that of a much larger pertur¬ 
bation (in Figures 14 and 15 ). This promotes the possibility that, in such a replication process, 
the timing of the error is equally or even more important than the size of the error. 


3.3 Mandelbrot sets for nonperiodic templates 

Due to the lack of regularity in the occurences of cq and ci in the template, the fractal dimension 
of the template Mandelbrot set seems to be lower than that for periodic templates. We conjecture, 
however, that the Hausdorff dimension is more uniform along the boundary of a set whose template 
exhibits less regularity - in contrast with the situation for regular templates. There, the fractal 
behavior was much more variable along the boundary of the corresponding Mandelbrot set, with 
portions with very high Hausdorff dimension, and portion with low dimension. 



Figure 10: Template Mandelbrot set for a random template, for fixed cq = 0. Each panel shows 
at higher resolution a zoomed in window, marked on the preceding panel. The stucture persists at 
higher and higher levels, suggesting the preservation of fractal features in the Mandelbrot set, even 
for random templates. 


4 Fixed-map Mandelbrod sets 

In this section, we suggest some alternative ways of viewing the fixed map Mandelbrot set, and 
of observing how the template affects the properties of the orbit Os(0), for a fixed parameter pair 
(co,ci) G C. As a book keeping method for truncated templates, we propose to consider each 
binary sequence of length L as equivalent to the binary representation of a real number 0 < n < 1 
with precision 2“^. In other words, we can consider, for each 0 < a < 1, its binary expansion up 
to its L-th binary digit, always choosing the infinite expansion (when it is the case). For fixed cq 
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and Cl, we can use each expansion as the symbolic template, and check whether the orbit Os( 0 ) 
is bounded or not. We can construct the function: F: [0,1] —)■ {0,1}, given by F{n) = 1, if Og is 
bounded, and F{n) = 0, if Os(0) is not bounded. In Figure 11 we show a few such representations, 
for various choices of the complex parameter pair (co,ci), connecting the discontinuous points in 
order to make the structure of the set more visible. It is clear that the behavior of the fixed map 
Mandelbrot set F“^(l) depends crucially on the parameter values. However, an analysis of its 
topological properties (accumulation points, density etc), or its measure in [ 0 , 1 ] remains open for 
further studies. 




1 


0 


Template n 


Figure 11; Fixed map Mandelbrot sets, for different parameter pairs (co,ci) G C. A. ci = 
-0.117 - 0.76i, CO = -0.5622 - 0.62i; B. ci = -0.5622 - 0.62i, cq = -0.5622 - 0.62i; C. ci = 
—0.117 — 0.76, Co = —0.75; D. ci = —0.75, co = —0.117 — 0.856i. We created binary expansions 
of lentgh L = 15 for the numbers in the unit interval [0,1], which we then used as the symbolic 
templates for the iteration process. 

Another interesting representation can be obtained as a hybrid of fixing the map and fixing the 
template. With a fixed cq, we can measure, for each different value of ci G C, how likely it is (i.e., 
for how many of all templates of a certain length) that the orbit Os(0) is bounded. In Figure 12 we 
show a few such slices in the ci complex plane, with the colors representing the number of templates 
for which Os(0) is bounded. 
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Figure 12: Hybrid Mandelbrot sets for different values of cq: A. cq = 0; B. cq = —0.75; C. 
Co = 0.375 + 0.333i; D. Cq = i. For each ci (represented in the complex plane), we used colors to 
illustrate how many templates of length L = 20 lead to a bounded orbit Os(0). The color spectrum 
goes from blue (low) to red (high). 

5 Discussion 

5.1 Comments and future work 

In this paper, we described a problem which is, to the best of our knowledge, a new extension 
of a well known traditional problem in discrete dynamics. We defined template iterations of two 
quadratic maps, and we used standard numerical algorithms to investigate the Julia and Mandelbrot 
sets of template systems. While this study is a first step in establishing the framework and phrasing 
some conjectures, a lot of theoretical work remains to be done, possibly exploring a brand new set 
of questions in discrete dynamics, which are likely to require new, different methods than the ones 
used in the traditional context of Fatou an Julia. 

The dynamic behavior in the case of one iterated map, and subsequently the connectivity type 
of classical Julia sets have been shown to be completely determined by the map’s critical orbits. 
While related approaches may still be appropriate to use in the case of alternating maps [ 6 ], and 
even in the case of higher template periods, this direction is not helpful when attempting to study 
the case of nonperiodic templates (which is in fact more plausible as a model for applications in 
the natural sciences). That is because the concepts of periodic or critical orbits are in this context 
no longer well defined, hence one has to start building a theoretical framework from its bare bones. 
In our current work we are exploring new methods, in particular expanding the question in a more 
abstract framework, obtained by interpolating continuously in the space of binary templates. 

In the meantime, in our computational work we are continuing to investigate extensions of 
the Julia set, in particular in the context of networks of interconnected complex quadratic maps. 
While pairs of coupled logistic maps have been studied before in both real and complex case, we are 
interested to study the coupled behavior in higher dimensional networks, with possible applications 
to understanding dynamics in neural networks. We are focused in particular on understanding 
the effects of the network architecture (e.g., graph Laplacian) and the dynamic properties of the 
ensemble as a whole, and on the topology of the “Julia set” of the networked system. 

5.2 Applications to genetics 

When a cell divides, it has to copy and transmit the exact same sequence of billion nucleotides to 
its daughter cells. While most DNA is typically copied with high fidelity (polymerase enzymes are 
amazingly precise when performing DNA synthesis), errors are a natural part of DNA replication, 
with rates of about 1 per 10^ (polymerases sometimes inserting too many or too few, or erroneous 
nucleotides into a sequence). Human diploid cells have 6 billion base pairs, and each cell division 
makes about 120,000 errors [ 12 ]. 
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Cells have evolved highly sophisticated DNA repair processes, aimed to promptly fix most of 
these errors. Some errors are corrected right away, during replication, through a repair process 
known as proofreading. Proofreading fixes about 99% of the errors, but that’s still not good 
enough for normal cell functioning. Some errors are corrected after replication, in a process called 
mismatch repair. Incorrectly paired nucleotides that still remain following mismatch repair become 
permanent mutations after the next cell division: once established, the cell no longer recognizes 
them as errors [12], passing them on to next generations of cells and (if the errors occur in gemetes) 
even to next generations of the organism. 

When the genes for the DNA repair enzymes themselves become mutated (the iterated function 
changes in the long-term), mistakes begin accumulating at a much higher rate. Mutation rates 
vary substantially among taxa, and even among different parts of the genome in a single organism. 
Scientists have reported mutation rates as low as 1 per 10®-10® nucleotides, mostly in bacteria, and 
as high as 1 per 10^-10^ nucleotides in humans [10]. Cells accumulate mutations as they divide. 
Even mutation rates as low as 10“^*^ can accumulate quickly over time, particularly in rapidly 
reproducing organisms like bacteria. 

Our system can be viewed as a theoretical framework for studying iterated replication mecha¬ 
nisms which are subject to errors at each iteration step. In genetics, polymerases replicate identi¬ 
cally the DNA strand at division, which in turn governs the development of the cell, and is passed 
on at the next cell division. In our model, the original system/cell (in our case, the complex z-plane) 
is programmed to evolve according to a certain sequence of steps, leading to emergence of some 
features and extinction of others. For example, an initial .^o which iterates to oo may represent 
a cell feature which becomes unsustainable after a number of divisions, while an initial .^o which 
is attracted to a simple periodic orbit may represent a feature which is too simple to be relevant 
or efficient for the cell. Then the points on the boundary between these two behaviors (i.e., the 
Julia set) may be viewed as the optimal features, allowing the cell to perform its complex function. 
An error at the level of the iteration function at one particular iteration step is equivalent to a 
mutation that occurred at one of the cell division steps. The new cell/complex plane is then used 
as template for the next iteration/division; one can study how the features of the cell are affected 
in the long term, when such an error passes undetected by the repair mechanisms. In our paper, 
we considered situations where such an occurrence is singular, random/occasional or periodic. It 
is clear that mutations accumulated over a long period of time may lead to serious changes in 
the structure of the later cells (different topological properties of the Julia set). Our model also 
addresses the timing when the errors occur, and illustrates how a mutation in the early iterative 
stages can lead to substantially more dramatic consequences on the result (Julia set) than the same 
error if it happens later in the process. 

The construction of mathematical models to help understand DNA replication and repair would 
be highly desired, since these are crucially important and complex mechanisms to study and under¬ 
stand. In eukariotic cells, accumulating mutations can lead to cancer. However, if DNA replication 
were perfect (mutation-free), there would be no genetic variation. Therefore, successful organisms 
had to construct optimal mechanisms, providing efficient DNA repair, but also enough variability 
for evolution to continue. A mathematical framework would be ideal for posing and contextualizing 
such questions. 
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Appendix 1: Effect of a propagating error along the binary template 





CO 



Figure 14: Effect of error propagation on the Julia set for the desired funetion with parameter ci = —0.117 — 0.856f (whose elassieal 
Julia set is shown in Figure le) and error parameter cq = 0 (whose elassieal Julia set is the unit eirele). The perturbation of to fc^ was 
introdueed sueeessively at the iterations A: = 1, 2, 3, 4, 5,10, 30 and 200 in a truneated template of length N = 200 (eaeh Julia set is represented 
in one of the figure panels, from left to right). 
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Appendix 2: Template Mandelbrot slices 



Figure 17: Mandelbrot slices for a lattice of cq or step 0.2, ranging from -0.6 to 0.6 in the real direetion, and from 0 to 0.8 in the 
imaginary direetion. The sliees are symmetrie with respeet to the real axis (not snown). 
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